Phonon-assisted tunneling regimes in diatomic molecules 
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Electronic transport in diatomic molecules (two-level systems) connected to metallic contacts is 
analyzed in the presence of competing electron-electron and electron-phonon interactions. We show 
that phonon emission and absorption processes are strongly modified when a Coulomb energy U is 
included, as the phonons open channels that can result in destructive or constructive interference 
effects. Resonance conditions for these processes produce dramatic effects both in the density of 
states at the molecular sites, as well as in the conductance through the system. We find in particular 
an enhanced Rabi-assisted tunneling due to phonons, as the resonance conditions are met, which is 
made more evident for increasing temperatures. These effects are controllable by voltage gating of 
the molecular sites, and should be accessible in current experiments. 

PACS numbers: 73.23.-b, 73.63.Kv, 71.38.-k, 73.63.-b 



I. INTRODUCTION 

Electronic transport properties of nanometer scale 
structures such as natural and artificial molecules (e.g. 
quantum dot structures) attached to metallic reservoirs 
have received a great deal of attention recently. Three 
main reasons drive the active research in the area: i) 
these systems have a fundamental role in the devel- 
opment of the field of "molecular electronics" where 
molecules serve as electronic components^ ii) they are 
relatively "simple" structures in which to address differ- 
ent conceptual issues regarding the nature of charge and 
spin transport) 2 - and iii) they represent the prototypical 
gateable structures at the nanoscale 4 

One of the most defining characteristics of these molec- 
ular systems is the strong spatial confinement that re- 
sults in enhanced electron-electron interaction (EEI), and 
profoundly modifies their transport properties. Well 
studied examples of the effect of EEI in confined sys- 
tems are Coulomb blockade and Kondo effect in quan- 
tum dots4»5iSi2iSiS Molecular systems may also exhibit 
strong electron-phonon interactions (EPI), which play 
a significant role in their physical properties. Several 
works have addressed the effects of EPI in quantum- 
dot molecules for different geometrical arrangements, in 
a perturbative regimei 10 ' 11 ' 12 It is known for example 
that phonons cause a broadening of Coulomb blockade 
peaks and the appearance of satellite peaks in the non- 
linear transport regime. Phonons also provide relaxation 
mechanisms (through inelastic scattering processes) that 
affect the electronic conductance, specially in molecular 
systemS ) 13 ' 14 i 15 as well as polaronic shifts of the electronic 
levels J£ 

While the effect of these different interactions on elec- 
tronic transport properties has been studied extensively 
in regimes where one of them is dominant, not much is 
known when both interactions become comparable. Re- 
cently, two works have addressed this issue. Al-Hassanich 
et al. have studied the effect of electron-phonon interac- 



tions in the transport properties of a single-level quan- 
tum dot connected to two metalic contacts in the Kondo 
regime^ They found an unexpected cancelation of the 
conductance resulting from interference between the pure 
electronic and phonon-assisted tunneling channels. 

Our group reported on the effect of competing interac- 
tions on the conductance of a diatomic molecule coupled 
to metallic reservoirs in the Coulomb blockade regime 
and for temperatures above the Kondo scaled That 
work focused on a regime where the energy difference be- 
tween the two local electronic levels is much larger than 
the broadening of the levels due to the coupling to the 
leads. In that case, indirect coupling between the local 
levels is negligible and interference of the different trans- 
port channels from one lead to another is not impor- 
tant. Utilizing the equations of motion for the various 
Green's functions involved, we found a splitting in the 
electronic density of states (DOS) of the molecule un- 
der suitable resonance conditions. This effect was iden- 
tified as Rabi-assisted tunneling, because of the formal 
similarity with the well-known Rabi splitting phenom- 
ena in atomic systems^ although in this case it is due 
to the interplay of EEI and EPI. The physical origin of 
the Rabi-assisted tunneling can be traced back to the ap- 
pearance of new channels for electron tunneling between 
the metallic reservoirs. These channels, due to EEI, pro- 
vide new resonance conditions for phonon emission and 
absorption. An immediate consequence of this result is 
the possibility of important changes in the conductance 
as function of the gate potential that controls electron 
occupation in the molecule, and which would be evident 
in experiments. This surprising result warrants a closer 
examination of the approximations used in the solution 
of the model and a thorough analysis of the various res- 
onance and interference conditions stemming from the 
EPI in the Hamiltonian. 

In the present paper, we address these issues and in- 
clude a complete analysis of the effects of the indirect 
coupling between the two electronic levels induced by the 
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FIG. 1: Schematic representation of the model system. Only 
the electron-phonon interaction connects the local sites (e° < 
ep) via phonons of frequency ljo and with coupling constant 
A. 



metallic leads. We now explore the regime where the lo- 
cal electronic level spacing tends to zero, which makes 
the phase effects much more important in the coherent 
transport processes. We also study in detail all resonant 
regimes and conducting channels provided by the EPI. To 
solve the equations of motion, we use a Hubbard-likc de- 
coupling procedure for all the Green's functions involved 
in the calculation.— This approach has proven reliable 
in the description of properties of interacting systems at 
temperatures above their characteristic Kondo tempera- 
ture. These calculations allow us to confirm the existence 
of the Rabi-splitting in the DOS and the appearance of 
resonances and anti-resonances not obtained previously. 
We show that there are structures in the DOS depend- 
ing upon resonance and anti-resonance conditions whose 
activation is level-occupation-dependent, reflecting the 
many-body nature of the problem. 

We further explore how the DOS structures and res- 
onance conditions affect the conductance through the 
diatomic molecule structure. We find that different 
Coulomb blockade peaks in the absence of EPI, appear 
split whenever phonons are present and this effect is en- 
hanced at higher temperatures. 



II. MODEL AND THEORETICAL METHOD 

Our model consists of a diatomic molecule with a single 
level per site (or equivalently a quantum dot pair) with 
local energies e° and A (we assume the non-degenerate 
case with e° > e°), as schematically shown in Fig.[TJ 

Each single-level site (or quantum dot) is connected 
independently to two external metallic leads. We in- 
clude local (on-site) electron-electron interactions and an 
electron-phonon inter-site interaction. As described the 
model can represent various experimental geometries; in 
particular, this model was studied in the limit of non- 
interacting electrons, to explore interference effects^ 
The total Hamiltonian is then written as 
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and the electron-phonon interaction is given by 
Hei-ph = <J a(7 cp a + h.c) , 



(3) 
(4) 

(5) 



where (b) creates (annihilates) a phonon with energy 
hoj , Ei = e° — eVg, and i = a, /?. We have assumed that 
U a = Up = U, and included independent gate voltages 
Vg that control the electron occupation of the molecular 
sites by shifting the localized energies with respect to the 
Fermi energy of the left and right electrodes. Also, 
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Hmoi-ieads connects the diatomic molecule to the leads; 
for simplicity, we consider t^ 7 = t' , for all 7 = R, L and 
i = a, (3. Away from the Kondo regime, the main effect 
of the leads is to broaden the energy levels of the dots 
through the tunnel couplings if . To describe the trans- 
port of electrons through the molecule, we calculate the 
local retarded Green's functions defined in the usual way 



as 
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where A and B are operators and (• • ■ ) denotes the ther- 
modynamic average (for T ^ 0) or the ground state ex- 
pectation value (for T = 0). In the frequency domain, 
the equation of motion for Green's functions is written 
as 



u,((A;B)) = ([A,B} ± ) + (([A > HUB)) ) 



(9) 



where the subindex in the brackets denotes commutation 
(— ) and anti-commutation (+) relations. The expression 
for the local electronic Green's function is: 



Gf» = ((c ia ;cl)), 



a, (3. 



(10) 



The assumption of identical leads allows the introduc- 
tion of symmetric and antisymmetric fcrmionic operators 
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After this transformation, H mo i and H m o\,-ie.ads read as 

Hi ea ds = ^ efcC 7fco .C 7 fecr. (12) 

-r=S,A 

and 

Hmol-leads = ^ f'c^CSfeo- + fl.C. , (13) 

where we have assumed epk = £Lk = £fe, and V = \^2t' . 
In the symmetric-antisymmetric basis, the coupling be- 
tween the molecule and the leads occurs only through 
the symmetric channel. The other terms in the Hamilto- 
nian are not affected by this transformation. We obtain 
the relevant Green's functions by using the equation of 
motion technique. As an illustrative example, let us de- 
scribe the steps taken in the calculation of the Green's 
function for level a : 

(uj - e Q )((c Qcr ;cJ lCT )) = 1 + U{{n aa c aa ]c\, a )) 
+A«&t C/Jo .;4 r » 

+ £?<(c sfcCT ;cL». (14) 

A; 

To obtain a closed set of equations, we apply a Hubbard- 
like decoupling procedure. In particular, for Green's 



In all these expressions g(e) refers to the Green's function 
of the leads at energy e (both leads have the same Green's 
function in equilibrium). In Eq. (|17|) . the expression for 

the self-energy E^(e) of an electron in level i, involves 
energies and occupation numbers of electrons in the other 
level j. These occupation numbers (nj a ) have to be cal- 
culated self-consistently. Since we are interested in the 
non-magnetic regime, the occupation numbers are taken 
to be spin-independent {(nj S ) — (nj a )) and they will be 
denoted just by (rij). Notice that the electron self-energy 
includes the indirect coupling between the two levels due 
to the presence of the leads through the second term in- 
side the brackets in the denominator of Eq. (fT7|) . The rel- 
evance of this term reduces as the inter-level separation 
in the dots increases. As a consequence, the net effect 
of the electron self-energy is to produce a finite width 



functions that involve operators for both levels a and 
(3 ("mixed Green's functions"), we close the equations 
by considering the following approximations: 

((tfJ aS c ps c aa -J a<J )) ps (6 t c^c (3er )((c a(T ;4 (7 ))(15a) 
((&4a c «sC Q a;cL)) ps (&4-c a ^}((c aCT ;4 CT )} (15b) 
((tfkwjcL)) ps (&t6>« CQCT ;4 CT )}. (15c) 

The first two equations are the standard Hubbard de- 
coupling procedure. At this level of the approxima- 
tion, correlations between electron operators with oppo- 
site spins are neglected. The third equation decouples 
the phonon mode from the electron operators, as sug- 
gested by Zubarev^ Finally, by keeping up to second 
order terms in A in the expression of the self-energies, we 
obtain: 

^(e) = { / £ " e ^;" £ 7^ 1 -^(e)}" 1 ! (16) 
y\e - a - U (1 - (n«r))] J 

with S iCT (e) = S^ ) (e)+S|f ) (e). Here S^ } (e) is the elec- 

tron self-energy in level i due to the leads, and ' (e) is 
the phonon-related self-energy. They have the following 
expressions: 



I 

in the DOS peaks, to shift their positions by an amount 
ReE^, and to give rise to transport interference effects 
due to the topology of the system. 

Equations (fTH|) and (fT9]l are the phonon self-energies 
calculated up to second order in the electron-phonon in- 
teraction strength A. They involve phonon occupation 
numbers (b^b), which for simplicity can be calculated 
from the free-boson Hamiltonian. Although this approx- 
imation is not essential, it is convenient to simplify the 
numerical calculation, and it does not introduce any ar- 
tifact in the results, as it is a frequency independent c- 
number. 

From Eq. (j 16[) . we obtain the total DOS by using 
p a {e) = (-l/7r)[lmG£ Q (e) + ImG^(e)]. In the limit 
where the electron-phonon coupling is negligible (A — > 0), 
the DOS presents four Coulomb blockade peaks coming 
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A 2 [(fctfr) + ( npa )] [e - e p + huj - U(l - (n ps ))} 

(e - ep + Hluq) (e - ep — U + huj ) - [e - ep + Hlu -U(l- (nps))] t' 2 g(e + hu ) ' 

A 2 [fib) + 1 - M] [e-e Q - huj Q - U(l - (n aS ))] 

(e - e a - huj ) (e - e a - U - tkj ) -[e-e a - huj Q -U(l- {n aS ))\ t' 2 g(e - hv ) 



(17) 
(18) 
(19) 
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FIG. 2: Schematic representation of the various resonances 
and the active levels. The figures correspond to the resonance 
conditions falBPal (b) and (c) [20b] and fdlEDcl 

from the original atomic levels located at e a , e a + U, 
ep and ep + U. As the electron-phonon interaction is 
turned on. these peaks show Rabi splittings, at char- 
acteristic energies, corresponding to phonon-absorption 
and/or phonon-emission processes at resonance condi- 
tions. These resonance conditions occur when the 
phonon energy Hljq matches the energy difference be- 
tween the two energy levels. These resonance conditions 
occur when 

fiw + U = Ae (20a) 
hoj Q = Ae (20b) 
tkj - U = Ae , (20c) 

where Ae = ep — e a . Figure [2] shows a schematic repre- 
sentation of the effective levels active in each resonance. 
Notice that the resonance conditions do not explicitly 
depend on the electron occupation number (unlike the 
simpler approximation results in (HI). However, when 
the system is found in one of these resonances the split- 
ting 8 a (x ReSi pft, ' ) for level a (and Sp oc ReSj 3 p ' i ' ) for level 
/3), depends on the gate voltage through the occupation 
numbers. 

Using the expressions above for the Green's functions, 
we proceed to calculate the effect of electron-phonon in- 
teractions on the DOS and the conductance of the sys- 
tem. Transport properties are calculated in linear re- 
sponse in the fully interacting regime (U ^ 0). The 
conductance has the following expression in the Keldysh 
formalism) 2 ^ 

G = 4 7 r 2 ^L(^)pH(^)|G I1 ( W )| 2 | w=eF - (21) 

where 

G u (w) = [§(w)tf [G aa + G P p + G a p + Gp a ] (22) 

is the Green's function that describes the dynamics of 
an electron with energy to from the site —1 to 1 (rep- 
resenting the leads, see Fig. [1} and g(u>) is the Green's 
function of the leads, calculated as disconnected from 



the molecule. As the four terms of Eq. (|22"|) are complex 
quantities with different phases, their contributions to 
the current can interfere constructively or destructively 
as we discuss below. The off-diagonal Green's functions 
(G a p and Gp a ) in the conductance are also calculated 
by equation of motion techniques at the same level of 
approximation described above. Finally, notice that we 
have neglected inelastic processes in the present calcula- 
tion because we assume that the system is essentially at 
thermodynamical equilibrium, at the small bias voltages 
of interest. 



III. NUMERICAL RESULTS 

Equations <P~^|) - (|19|) are solved numerically and used to 
calculate the DOS and conductance in (|2"Tj) , as described 
above. For convenience, the Fermi energy is set to zero 
(cf = 0). In what follows, all the physical quantities are 
given in units of the inter-level separation Ae. The levels 
e° and ejg are set at 1.0 and 2.0, respectively (above ep)- 
The effect of an increase in the gate voltage V g is lowering 
these levels (initially towards the Fermi energy). In the 
next few figures, we assume equal gating of both sites, so 
that V g a = V^ = V g . 

Let us consider first the results for the interacting sys- 
tem in the absence of electron-phonon ineraction (U ^= 
0;A = 0). Figure © shows the results for the DOS 
(top) and the conductance (bottom) for the particular 
values U = 0.4 and T = 0.0025 (energies in units of 
Ae = ep — e Q ). As expected, the conductance exhibits 
peaks located at gate voltages near e Q , e a + U, ep and 
ep + U; these are the well-known Coulomb blockade peaks 
that appear whenever each level is in resonance with the 
Fermi energy. Notice slight shifts from the anticipated 
energy values, caused by the non-vanishing real part of 
the electron self-energy discussed above. 

Figure 2] shows the DOS as a function of gate voltage 
V g and energy w in the presence of electron-phonon in- 
teraction. The different plots refer to different values of 
the phonon energy Hljq. In all cases, the electron-electron 
interaction is set at U — 0.4, the electron-phonon cou- 
pling strength A = 0.2, and T = 0.0025. The values for 
the phonon energies, frio® = 0.6,1.0,1.4, have been cho- 
sen to fit the resonance conditions given by Eqs. (|20aD , 
(|20b[) and (|20c|) . respectively. From the figure, it can 
be observed that in the first resonance condition for the 
smaller tkoo (= 0.6, top panel) the peaks correspond- 
ing to the levels e a + U and ep are split, mixed by the 
phonon channels. The splitting for the peaks correspond- 
ing to the level e a + U is w 28 a (« at a gate voltage 
V g w 1.4), while the splitting for the level ep is w 25 p 
(« 0.05 at V g as 2.0, for example). Notice that an analy- 
sis of Eqs. (fT7|) -([T^ l) . reveals that when level e Q coincides 
with eF the splitting ks 2S a becomes small. This can be 
understood in terms of the occupation number for the 
level P which in these conditions remains close to zero 
(this is seen by analyzing the value of the ReS Q , which 
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FIG. 3: (Color online) Density of states map as function of 
gate voltage and energy (top panel) and conductance as func- 
tion of gate voltage (bottom) for a molecule with electron- 
electron interaction, but no electron-phonon interaction. In 
units of Ae = ep — e a , the parameters are A = 0, U = 0.4, 
Hivo = 0.6, f = V2/5, e° = 1, and T = 0.0025. Notice energy 
scale in top panel is shifted by e a . 

approaches zero) . A similar analysis shows that when the 
level ep crosses ep the splitting 26 p is also small, since 
Re£^ sa 0, a consequence of the levels e a and e a + U 
being far below ep (almost full levels). The splitting of 
the levels in the DOS is clearly, and most importantly, 
manifested in the conductance as a function of the gate 
voltage V g , as shown in Fig. [5] 

In the case of the second resonance condition, hujQ = 
1.0 (middle panel in Fig. 2J), there are two pairs of levels 
involved: e a , ep and e a + U, ep + U. The effect of the res- 
onance conditions on the conductance results to be much 
less important than in the previous case, as the splitting 
in the DOS at the Fermi level is much smaller here. No- 
tice, however, that in this case all four conductance peaks 
exhibit a certain degree of splitting (see middle panel in 



Fig. [5J), indicating that the interference effects involve 
two pairs of levels (as in Fig. [2b andHfc). 

For the last resonance condition, hwo — 1.4, the levels 
involved are e a and ep + U . These levels are quite far 
apart, separated by Ae + U, which results in the split- 
tings in the DOS and the conductance to be barely ob- 
servable, except for the sharp dip/splitting in the first 
and last Coulomb blockade peaks (lower panel, Fig. [SJ. 
We notice that the Rabi splitings are generally less im- 
portant for the last two resonance conditions. The reason 
for this behavior is that larger hujo means that the orig- 
inal levels participating in the phonon processes are far 
apart, decreasing the spliting energies 6 a and 6p. This 
is clearly reflected in the respective conductance curves, 
where the peak splittings are very small. 

We have also analyzed the dependence of these results 
on temperature. The top panel in Fig. [5] shows the con- 
ductance as a function of gate voltage for three different 
temperatures, T = 0.0025, 0.025, and 0.1. A relatively 
weak dependence on temperature is observed, consistent 
with the fact that the resonance conditions are not ex- 
plicitly dependent of the electron occupations. However, 
the self-consistency of the calculation induces an indirect 
occupation dependence, resulting in splittings that be- 
come more pronounced for increasing temperature (i.e., 
larger phonon numbers). 

In order to fully probe all the resonance conditions, we 
fix hujQ and U and control the energy spacing Ae by gat- 
ing only one of the localized levels. Figure [6] shows a color 
map of the DOS as function of frequency and energy of 
the variable level ep. In this plot ep is shifted down from 
above to below ep as V g increases; e a is kept fixed at a 
point just below the Fermi energy, e Q = —0.001, so that 
A e also changes with V g . Since the energy spacing is not 
constant, we choose instead A = D/8 as the energy unit, 
where D is the bandwidth in the leads. Let us analyze in 
detail the various peaks in this figure. Following a ver- 
tical line, such as ep = 2.25A, one finds DOS peaks at 
£a + fo^o an d e a + U (sa 0, 1, and 1.4), respectively. 
There is one more peak at e a + U + Hluq (sa 2.4), which 
is out of range of the vertical axis in this panel. Note 
that peaks at both e Q sa and e a + U appear since the 
site a is in the mixed valence regime. Now following a 
horizontal line, at u> = —A, we find peaks corresponding 
to the levels ep — hujQ, ep, ep + U— thu^ and ep + U, respec- 
tively One notices anticrossings of levels whenever one 
of the levels is a phonon replica, and crossings when both 
levels are "purely" electronic. The anticrossings reflect a 
direct interaction between the electronic levels mediated 
by phonons. The crossings appear because the levels in- 
teract only indirectly through the leads, an interaction 
that is energy dependent. 

The bottom of Fig.[6]shows the conductance of the sys- 
tem. Note the appearance of antiresonances in the con- 
ductance whenever a crossing or anticrossing coincides 
with the Fermi level. The antiresonances localized near 
1 and —0.4 are due to the anticrossings in the density of 
states discussed above. They result from the absence of 
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FIG. 4: (Color online) Map of the DOS as function of gate 
voltage and energy in the EEI and EPI competing regime. 
From top (in units of Ae), with fttJo = 0.6, huio = 1 and 
huj Q = 1.4. In all panels, U = 0.4, t' = y/2/5, A = 0.2, e° = 1, 
and T = 0.0025. The energy axes are shifted by e a , so that 
level a appears always at zero. 



level e Q , as it is always kept at the Fermi level. The full 
uninhibited a-level conductance is evident for ep > 1.5, 
where the /3-level is empty. Notice yet, that the total 
conductance exhibits dips and even complete vanishing 
due to the complex pattern of interference arising from 
competing channels of electrons traversing the molecule 
through both arms of the structure whenever ep is near 
the Fermi level. In particular, the vanishing conductance 
at ep ~ —0.5 appears from the strong level anticrossing 
lying right at the Fermi level (see upper panel in Fig. [5]). 

In Fig. [7] we plot the conductance of the system as 
function of ep, where the level a is shifted to a posi- 
tion below the Fermi energy (e Q = —0.1 A), in such a 
way that it essentially does not contribute to the con- 
ductance. As we apply a gate voltage to change ep, we 
observe a conductance curve consisting of three peaks. 
Two of them, located near ep = and ep = — 1.4A are 
the usual Coulomb blockade peaks, which are essentially 
not affected by the EPI. The effect of phonons, how- 
ever, can be seen most clearly in the appearance of an 
additional peak near — 0.5A, due to phonon-assisted pro- 
cesses. This peak arises from the anticrossing of the levels 
e a and ep + U — tkJo, opening a new channel for electron 
transport. 



IV. CONCLUSION 

In summary, we have studied in detail a diatomic 
model when both EEI and EPI are taken into account. 
We have presented an accurate description of the Rabi- 
assistcd tunneling phenomena studying the Green's func- 
tions of the system in the Coulomb blockade regime above 
the Kondo temperature. We have obtained well defined 
resonance conditions for the system. A detailed analy- 
sis of the various resonance conditions and interacting 
regimes shows that the most dramatic changes occur for 
weak EEI (U, hojo < Ae), since the EPI is enhanced when 
the states involved in the phonon emission and absorp- 
tion are close in energy. These results emphasize the rich 
physics involved in phonon-assisted phenomena, which 
we believe should be accssible in experiments. The fact 
that the phonon-induced resonance conditions are sub- 
tly temperature dependent makes the effect strong and 
sensitive to be probed in experiments. 
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FIG. 5: (Color online) Conductance as function of gate volt- 
age for system as in Fig. [4] From top to bottom, huio = 0.6, 
Tiujo = 1 and huo = 1.4; parameters as above. Top panel 
also shows T dependence of conductance for /two = 0.6; no- 
tice splittings in middle two peaks increase with increasing 
temperature. In the top panel the curves are offset for clarity. 
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FIG. 6: (Color online) Top: Map of the density of states as 
function of energy (vertical axis) and level position ep (hor- 
izontal axis). Bottom: Conductance as function of energy 
position ep. Energies in units of A = D/S. Other parame- 
ters are Huo = 1 ■ A, and U = 1.4A, i' = V2/5, A = 0.2, 
and T = 0.0025A. Level a is kept fixed at e a = -0.001A, 
essentially at the Fermi level and always contributing to the 
conductance and interference. Notice scale runs opposite 
each other in both panels. 
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FIG. 7: Conductance as function of position of level ep for the 
same case as in Fig. [SJ but with the e a level well below the 
Fermi energy, e a = —0.1/ Delta, so that it does not contribute 
to the conductance. 



